Breakdown of Energy Equipartition in a 2D Binary Vibrated Granular Gas 
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We report experiments on the equipartition of kinetic energy between grains made of two different 
materials in a mixture of grains vibrated in 2 dimensions. In general, the two types of grains do 
not attain the same granular temperature, Tg = l/2m < v'^ >. However, the ratio of the two 
temperatures is constant in the bulk of the system and independent of the vibration velocity. The 
ratio depends strongly on the ratio of mass densities of the grains, but is not sensitive to the 
inelasticity of grains. Also, this ratio is insensitive to compositional variables of the mixture such 
as the number fraction of each component and the total number density. We conclude that a single 
granular temperature, as traditionally defined, does not characterize a multi-component mixture. 

PACS numbers: 81.05.Rm, 83.10.Pp, 05.20.Jj 



When two different gases are put in thermal contact 
with each other they eventually reach thermal equilib- 
rium, with equipartition of kinetic energy between the 
gas molecules and a single temperature for the system as 
mandated by the zeroth law of thermodynamics. In this 
article we present an experiment designed to test the va- 
lidity of energy equipartition in a mixture of two kinds of 
macroscopic grains driven into a non-equilibrium steady 
state by external vibration. Kinetic theories for granular 
systems have long used the idea of a granular tempera- 
ture defined as the average kinetic energy of the particles, 
Tg — l/2m < >; the value of Tg in a given system is 
determined by a balance between the power input by ex- 
ternal driving sources and the energy lost in interparticle 
collisions. The results of such calculations of transport 
in single-component systems have been validated in both 
simulations ||l|, ||] and experiments It is natural to 
ask whether the variable Tg obeys a zeroth law, and in 
particular, whether a mixture of grains is characterized 
by a single, shared, granular temperature. 

The earliest theoretical developments of kinetic theo- 
ries for binary mixtures of spheres ^ or of discs 
start with a clearly identified assumption that a single 
temperature variable, Tg, characterizes the entire mix- 
ture. Most subsequent work preserves this assumption. 
An important exception is the work of Garzo and Dufty 
who study theoretically the "cooling" of a binary 
granular mixture from an initial distribution of parti- 
cle velocities. They conclude that the cooling rate of 
each species in the mixture is the same, although their 
temperatures are different throughout the cooling pro- 
cess. Experiments and simulations on binary mixtures 
of granular materials have primarily focused on segrega- 
tion and mixing of different types of grains j8|, and for 
the most part, have not considered the issue of equipar- 
tition. However, an experiment on a vibrated monolayer 
of grains by Losert et al. M found that grains of two 



different sizes and masses do not attain the same value 
of Tg. In this article, we show that such differences per- 
sist in the bulk of a 2D system and that in general, the 
two components of a mixture arrive at a ratio of granu- 
lar temperatures, 7 7^ 1. We explore the dependence of 
the temperature ratio, 7, on particle properties as well as 
on compositional parameters of the mixture. Our data 
argue for a redefinition of the granular temperature in 
order to accommodate extension to a multi-component 
system. 

We make binary mixtures of spherical glass balls with 
aluminum, steel or brass balls. The two relevant material 
properties being varied are mass and inelasticity (see Ta- 
ble |), and the two compositional parameters being var- 
ied are the number fraction of each component, x, and 
the number density, specified in terms of the the average 
occupied area fraction, pavg, of the particles in the cell. 
To determine the importance of particle material prop- 
erties we hold average number density and number frac- 
tion fixed and compare mixtures of steel (heavy-elastic) 
and glass (Hght-elastic), brass (heavy- inelastic) and glass 
(light-elastic), and aluminum (light-inelastic) and glass 
(light-elastic) in terms of the Tg of each component and 
the temperature ratio 7. Next, to determine whether fail- 
ure of equipartition is a bulk effect, we hold the number 
fraction fixed for a steel-glass mixture and measure the 



TABLE I: Some material properties of the spheres used in 
the experiment. 



Particle 


Mass 


Efi^ective 


Mass ratio 




[mg] 


inelasticity " 


w/ glass 


Glass 


5.24 


0.17 




Alum 


5.80 


0.31 


0.92 


Steel 


15.80 


0.21 


0.33 


Brass 


18.00 


0.39 


0.28 
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"The effective inelasticity quoted is the fraction of energy lost in 
a collision in single component system, averaged over all collisions. 
Since relative velocity of the spheres, the impact parameter, and 
rotational motion all affect the inelasticity of a given collision, the 
number quoted above is the mean of a wide distribution. 
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FIG. 1: (a) Density profile, p, of alum-glass {pavg ~ 0.096) 
versus height, z (in d). (b) Density profile, p, of steel-glass 
{pavg = 0.096), and brass-glass {pavg = 0.049, scaled) versus 
height, z. Dashed lines represent the limits of the vibrating 
walls, r = 56 g and x = 1/2 for all 3 pairs. 



temperature ratio as the number density of balls is var- 
ied. Finally, to test the dependence of the temperature 
ratio on the relative fraction concentration, we hold the 
number density fixed for a brass-glass pair, and measured 
the temperature ratio of the components as the number 
fraction of brass is varied. 

The balls (diameter: d — 1.600 ± 0.002 mm) are con- 
fined to a vertical, rectangular, Dclrin cage (32 d high x 
48 d wide x 1.1 c? thick) sandwiched between two parallel 
plates of anti-static coated Plexiglas. An electromechan- 
ical shaker (LDS 500L) vibrates the cage vertically at 
a frequency of 60 Hz and amplitudes up to 2.4 d, pro- 
ducing maximum accelerations, F, and velocities, uq, of 
56 g and 1.45 m/s respectively. The motions of the balls 
are recorded with a monochrome high-speed camera (Ko- 
dak Motioncorder) at a rate of 2000 frames/s. Each ball 
is located with a precision of 0.03 d (~ 40 iim). Glass 
balls are unambiguously distinguished from their metallic 
partners by using a combination of reflected and trans- 
mitted light, and by thresholding appropriately. Most 
of the results we discuss here are taken in a rectangular 
(10 d x 20 d) window centered in the geometrical center 
of the cell. 

The data presented in this article are all taken at 
relatively large values of excitation, corresponding to 
F > 32 g and vq > 0.86 m/s. Beyond this scale of excita- 
tion, the vertical density profile does not vary much and 
tends to an asymptotic profile that is symmetric about 
the center of the cell, as shown in Fig. ^ The density 
profile is identical for two different species when their 
masses are matched, even if they have different inelas- 
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FIG. 2: (a) Average granular temperature < Tg > versus 
squared vibration velocity Vq for alum-glass, brass-glass, and 
steel-glass mixtures. Open symbols and filled symbols repre- 
sent Tg for glass and the metal partner, respectively, pavg = 
0.096 and x = 1/2 for all mixtures, (b) Ratio between glass 
temperature and partner temperature, 7 = Tgiass/Tpartner, 
versus the squared vibration velocity, vq- 



ticity, as in the case of aluminum-glass mixture shown 
in Fig. |l|(a). In contrast, the profiles of the other two 
mixtures (steel-glass and brass-glass) in Fig. |l|(b), show 
a slightly higher concentration of heavy particles in the 
center of the cell and a more uniform distribution of light 
particles throughout the cell. However, in all three mix- 
tures, while there are differences in the overall density 
distribution, the particles in the mixture are locally well- 
mixed, with no obvious tendency for balls of a particular 
type to cluster |11 when vigorously excited as in this 
case. Finally, we have not observed any gradients of ve- 
locity or density in the horizontal direction. 

The principal observation of our experiment is that the 
two components of a binary mixture do not always equi- 
librate to the same granular temperature. Figure |^(a) 
shows Tg of three pairs of materials as a function of the 



squared vibration velocity of the cell, Vq, with the num- 
ber density being held fixed at pavg = 0.096 and the 
number fraction of glass balls fixed at x = 1/2. The 
granular temperature Tg is proportional to the squared 
vibration velocity Vq with the proportionality constant 
being a function of the inelasticity of the balls and num- 
ber density Pavg- As shown in Fig. |^(b), Tg equilibrates 
for the alum-glass mixture, but does not equilibrate for 
the steel-glass and brass-glass mixtures. In terms of 
particle properties, equipartition is achieved when the 
mass ratio is close to one (mgiass /malum — 0.92), and 
fails when the masses are mismatched {mgiass/^steei = 
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FIG. 3: Vertical temperature profiles for a steel-glass mix- 
ture. Dotted and dot-dash lines represent Tsteei and Tgiass 
and the full line represents the ratio between the glass and 
steel temperatures. Dashed vertical lines represent the limits 
of the vibrating walls. For this system pavg = 0.096, x = 1/2 
and r = 56 gr. 



0.29, TTT-giass / fTT-brass — 0.33). Even wheii the temper- 
atures of the two components do not equilibrate, the 
ratio of the temperatures do not vary as a function of 
Vq within the error bars, as shown in Fig. §(b). Fur- 
thermore, the ratio of temperatures are very close to 
one another [jg^s = Tgiass/Tsteei = 0.66 ± 0.06 and 
7g_h = Tgiass/Tbrass = 0.69 ± 0.09), cvcu though the 
inelasticity of the metal partner in each case is quite dif- 
ferent (Table |). We conclude that the temperature ratio 
differs from one when the mass ratio departs significantly 
from one, but does not depend sensitively on the relative 
value of the inelasticity of the two components. 

The failure of equipartition is not merely a boundary 
effect, nor is the value of 7 that we report in Fig. ^(b), an 
inhomogeneous average over different regions of the cell. 
This is shown in Fig. |^ where the temperature profiles 
of steel and glass are plotted as a function of the height 
of the cell, z[d], for vq = 1.45 m/s. In a collision with a 
moving wall which is effectively infinitely massive, both 
species pick up the same velocity, which implies that the 
heavier species (steel) picks up greater kinetic energy. 
The temperature ratio 7 is thus different from one at the 
walls; this is analogous to the previous observations of 
Losert et al [||. Remarkably, however, the temperature 
ratio quickly reaches a constant value (7 = 0.69) in the 
interior of the cell, even as the Tg of each component 
continues to decrease toward the middle of the cell as a 
result of losses from inelastic collisions. Therefore the 
effect of inelastic collisions is to drive the value of 7 in 
the interior of the cell toward a value different from one, 
which as we mentioned above, depends on the mass ratio 
of components, but is insensitive to their inelasticity, and 
to the value of vq. 

We pursue further the possibility that the value of 7 
obtained is due to boundary effects by increasing the 
number density, Pavg, of a steel-glass mixture holding 
the number fraction fixed at x = 1/2. While this is not 
exactly equivalent to increasing the system size, it does 



tune the ratio of system size to mean free path, so that 
we vary the number of collisions experienced by balls 
in transferring energy from the boundary to the inte- 
rior. The temperature ratio, 7 = Tgiass/Tsteei is plotted 
against of the square vibration velocity, Wq, in Fig. ^, for 
mixtures whose number density, pavg, varies by a factor 
of 4. Once again, the ratio of the two temperatures are 
essentially independent of Vq. Importantly, the tempera- 
ture ratio maintains an average value of 7 = 0.60 ± 0.06, 
with no obvious systematic trend as pavg changes. This 
result further confirms that the failure of equipartition 
is a bulk effect, and also suggests that the value of 7 is 
insensitive to the number density of the mixture. 

Finally, we study the effect on equipartition of vary- 
ing number fraction in a mixture at fixed density. Figure 
H shows the temperature ratio, 7 = Tgiass/Tbrass, of a 
brass-glass mixture as a function of Wg, for three differ- 
ent number fractions of brass, x: 1/2, 1/4, and 1/8. The 
temperature ratios for the different relative concentra- 
tions are all the same within experimental errors. This 
demonstrates that the components do not reach thermal 
equilibrium, even though the relative number of colli- 
sions with glass balls increase for a given brass ball, as 
the number fraction of brass is decreased. It has been 
predicted ||l^ that equipartition does not generally oc- 
cur even in the tracer limit of a single foreign particle in 
a sea of other particles. This is consistent with our ob- 
servation that the temperature ratio 7 does not depend 
on number fraction. 

As mentioned earlier, Garzo and Dufty study a 
freely-cooling binary system, and find a temperature ra- 
tio dependence on material parameters as well as on com- 
positional parameters. Since there are no comparable 
predictions for a driven system, we have compared our 
data with these predictions and find quite good agree- 
ment with our results as a function of material parame- 
ters. However, they predict a strong dependence on num- 
ber fraction of the components, where none is observed. 
The comparison between theory and experiment encour- 
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FIG. 4: Temperature ratio, 7 = Tgiass/Tsteei, in a steel- 
glass mixture with number fraction x = 1/2, plotted against 
squared vibration velocity, Vq. Different markers represent 
different number densities of the mixture. 
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FIG. 5: Ratio between glass temperature and brass tem- 
perature, 7 — Tgiass/Tbrass, versus squared vibration velocity 
of the cell, Uq. Different markers represent different num- 
ber fractions of brass for the same total number of particles 

{Pavg = 0.049). 



nents does robustly satisfy some general trends. The 
temperature ratio is constant in the interior of a system, 
and is insensitive to vibration velocity and inelasticity, 
while depending strongly on mass density. Furthermore, 
the ratio is insensitive to compositional parameters such 
as number fraction of each component and average num- 
ber density. This opens the possibility that just as in 
the work of Garzo and Dufty |7| for a cooling system, 
there is a well-defined and calculable ratio of tempera- 
tures achieved in a steady state system. Thus a modified 
definition of granular temperature with a prefactor de- 
pending chiefly on material properties, may adequately 
describe a multi-component system. Numerical simula- 
tions could play a useful role, in that material properties 
may be continuously tuned, as opposed to experiments in 
which only a limited pool of materials is readily available. 
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In conclusion, we have observed that energy equiparti- 
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lar system. However, the ratio between the steady state 
values of the granular temperatures of the two compo- 
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